The Angular Clustering of Galaxy Pairs 
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ABSTRACT 

We identify close pairs of galaxies from 278 deg^ of Sloan Digital Sky Survey commissioning 
imaging data. The pairs are drawn from a sample of 330,041 galaxies with 18 < r* < 20. We 
determine the angular correlation function of galaxy pairs, and find it to be stronger than the 
correlation function of single galaxies by a factor of 2.9 ± 0.4. The two correlation functions have 
the same logarithmic slope of 0.77. We invert Limber's equation to estimate the three-dimensional 
correlation functions; we find clustering lengths of tq — 4.2 ± OAh^^ Mpc for galaxies and 
7.8 ± 0.7 h^^ Mpc for galaxy pairs. These results agree well with the global richness dependence 
of the correlation functions of galaxy systems. 

Subject headings: galaxies: evolution - galaxiesiformation - galaxiesdnteractions clusters:groups 
1. Introduction 
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Clusters are observed to be more strongly clus- 
tered than galaxies (Bahcall & Soneira 1983; 
Bahcall 1988; Postman et al. 1992; Croft et al. 
1997; Abadi et al. 1998). Moreover, the clus- 
tering strength of groups and clusters has been 
shown to increase with richness (Bahcall & Soneira 
1983, Szalay & Schramm 1985, Bahcall & West 
1992). Groups of galaxies have a smaller correla- 
tion length than that of rich clusters (Merchan et 
al. 2000) but larger than that of individual galaxies 
(Connolly et al. 1998, Loveday et al. 1996, Infante 
& Pritchet 1995). A number of authors have ex- 
amined the correlation function of groups of galax- 
ies from relatively shallow redshift surveys (Jing 
& Zhang 1988; Maia & da Costa 1990; Ramella 
et al. 1990; Trasarti-Battistoni et al. 1997; Gi- 
rardi, Boschin & da Costa 2000; Merchan, Maia, 
& Lambas 2000), finding that the groups exhibit a 
correlation function somewhat stronger than that 
of galaxies. The richness dependence of the corre- 
lation function is generally explained in terms of 
high-density peak biasing of the galaxy systems 
(Kaiser 1984), and is seen in cosmological simu- 
lations {e.g. , Bahcall & Cen 1992, Colberg et al. 
2000 and references therein). 



1 



This paper quantifies the angular clustering of 
pairs of galaxies, thus exploring the dependence 
of clustering on system richness in the regime be- 
tween single galaxies and groups. In future work, 
as an appropriately large sample becomes avail- 
able, we will investigate the clustering of groups 
of galaxies (see also Lee & Tucker 2001). 

The Sloan Digital Sky Survey (SDSS; York et 
al. 2000) commissioning data (Stoughton et al. 
2001) provide a photometrically reliable catalog of 
galaxies over a large field to r* « 22.0. We define 
a uniform sample of compact galaxy pairs from 
these data over an area of 278 dcg^. Scranton et 
al. (2001) show that systematic effects on galaxy 
clustering due to star-galaxy separation, varying 
seeing, photometric calibration, reddening, and so 
on, are small in these data. The resulting galaxy 
clustering analysis is presented by Connolly et al. 
(2001), Tegmark et al. (2001), Szalay et al. (2001), 
and Dodelson et al. (2001). In the present pa- 
per, we extend these results by comparing the an- 
gular correlation function of galaxies to that of 
galaxy pairs. Our magnitude slice samples galax- 
ies with a median redshift of roughly 0.22, appre- 
ciably deeper than that of previous group studies. 

Assuming an appropriate redshift distribution, 
we invert Limber's equation to determine correla- 
tion lengths for both galaxies and pairs, and its 
dependence on the number density of these sys- 
tems. 

In § 2 we describe the Sloan Digital Sky Survey 

imaging data, the properties of the galaxy cata- 
log, and the definition of galaxy pairs. Estimation 
of the angular correlation function of galaxies and 
of pairs is presented in § 3. In § 4 wc invert the 
angular function to determine the spatial correla- 
tion scales of galaxies and pairs. The dependence 
of the correlation scale on richness is presented in 
§ 5, and we give our conclusions in § 6. 

2. The Sloan Digital Sky Survey Imaging 

Survey 

The Sloan Digital Sky Survey is a photometric 
and spectroscopic survey of 1/4 of the sky, above 
Galactic latitude of ~ 30° (York et al. 2000). The 
photometric data are taken with a dedicated 2.5 m 
altitude-azimuth telescope at Apache Point, New 
Mexico, with a 2.5° wide distortion-free field and 
an imaging camera consisting of a mosaic of 30 



imaging 2048x2048 SITc CCDs with 0.4" pixels 
(Gunn et al. 1998). The CCDs are arranged in six 
columns of five CCDs each, using five broad-band 
filters (u, g, r, i and z). The total integration time 
per filter is 54.1 seconds. Each column of CCDs 
observes a scanline on the sky roughly 13' wide; 
the six scanlincs of a given observation make up a 
strip, and two interleaved strips give a stripe 2.5° 
wide. The measured survey depth at which repeat 
scans show 95% reproducibility is 22.0, 22.2, 22.2, 
21.3, and 20.5 magnitudes for the 5 filters, respec- 
tively. The SDSS photometric system is measured 
in the AB^ system (Oke & Gunn 1983, Fukugita 
et al. 1996). 

Several aspects of the SDSS photometric pipeline 
(Lupton et al. 2001) are worthy of mention in the 
context of studies of pairs and groups of galaxies. 
We use Petrosian (1976) magnitudes, as described 
by Blanton et al. (2001), Yasuda et al. (2001), 
and Stoughton et al. (2001). Overlapping images 
are deblended consistently in five colors, using 
an algorithm which makes no assumption regard- 
ing profile shape or symmetry thereof, and which 
works with an arbitrary number of overlapping 
objects. Visual inspection shows that the mor- 
phologies and photometry of overlapping galaxies 
are correctly determined for pairs of objects of 
similar brightness whose centers are separated by 
as little as 3 arcsec. 

Star-galaxy separation is accurate at the 99% 
level to at least r* = 20.51"*, and better than 90% 
for r* < 22 (for < 1.5" seeing); we discuss this 
further below. The data are deep enough that the 
typical photometric error in the Petrosian magni- 
tude of an r* = 20.8 galaxy is 0.05 mag, plus 0.03 
magnitude error to be added in quadrature due to 
uncertainties in the overall photometric calibra- 
tion. 

2.1. The Data 

We use imaging data taken during the commis- 
sioning period of the SDSS (on 21 and 22 March 
1999), which together make up a stripe 2.5° wide 
and 100° long centered on the Celestial Equator 
(runs 752 and 756); these data are included in 
the SDSS Early Data Release (Stoughton et al. 



Because the photometric cahbration is not finalized, we 
refer to observed photometry with asterisks. See Stoughton 
et al. (2001) for a full discussion. 
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2001). These runs lie within the area 1^.1 < 
R.A. (2000) < 16''.8 and -1°.26 < Dec. (2000) < 
1°.26, although the two strips overlap for only the 
central part of this right ascension range. The see- 
ing ranged from 1.2" to 2.5". 

For each of the six scanlines of the two runs, we 
use the data that is more than 30" from the edge 
of the CCDs. We carried out tests of the robust- 
ness of the derived angular correlation function to 
the masking of bright stars, and found it to be in- 
sensitive to the masking; in what follows, we do 
not mask the stars. We carry out the pair counts 
on each of the twelve scanlines separately, normal- 
izing in each one, and then average the results. 
Note that this differs from the approach of Scran- 
ton et al. (2001), who analyze the twelve scanlines 
together; they also define a mask which excludes 
bright stars and regions of particularly poor see- 
ing. We will sec below that our angular correlation 
function is essentially identical to theirs. The re- 
sulting sample has an effective area of 278.13 deg^. 

The turnover in the galaxy number counts as 
a function of magnitude (Yasuda et al. 2001) pro- 
vides a good estimate of completeness. The galaxy 
number comits in r* rise as dlogN/dm w 0.4 up 
to i? < 25 (Infante, Pritchet, & Quintana 1986; 
Tyson 1988). Thus any turnover relative to this 
power law at r* < 25 can be attributed to in- 
completeness. Figure 1 shows the galaxy counts 
for our sample, after correcting for reddening us- 
ing the extinction maps of Schlegel, Finkbeiner, & 
Davis (1998) (sec Yasuda et al. 2001 for a thor- 
ough discussion of the SDSS galaxy coimts). The 
thin solid lines represent the counts for each of 
the 12 scanlines analyzed. The solid dots arc the 
overall counts for the sample. The bold solid line 
is the best power law fit to the data points. The 
slope of the counts is dlogN/dm — 0.46. The 
counts begin to turn over noticeably fainter than 
r* = 20.5 (although Scranton et al. 2001 show 
that one can push star-galaxy separation fainter 
for galaxy clustering analyses). We therefore take 
our galaxy catalog to be complete to r* < 20.5. 

2.2. The Catalogs 

The next step is to generate catalogs of galaxies 
and galaxy pairs. 

We wish to work in a fairly narrow magnitude 
range to minimize redshift projection effects. We 



also want to avoid problems with completeness, 
and work in a range where the redshift distri- 
bution of galaxies is well- understood (see below). 
With this in mind, we define a sample in the range 
18 < r* < 20, after correction for reddening as 
above; this includes 330, 041 galaxies over an area 
of 278.13 deg^. The distribution on the sky of the 
galaxies in our sample is shown in Figure 2, which 
shows the region centered on the right ascension 
range in which the two strips overlap. Although 
these data are drawn from 12 separate scanlines, 
the scanline boundaries are invisible where they 
overlap. 

The reliability of the catalog depends on pho- 
tometric uniformity, the performance of the ob- 
ject deblender, the number of spurious detections 
around bright stars and edge effects. These issues 
are discussed by Scranton et al. (2001) and Lup- 
ton et al. (2001). In brief, overlaps between ad- 
jacent scanlines show photometric zero-point off- 
sets of less than 0.03 mag; similar conclusions are 
reached from the agreement in galaxy counts be- 
tween the scanlines as shown in Figure 1. Vi- 
sual examination of pairs and triplets of galaxies 
found in the sample (as described below) shows 
that the deblender works well in over 95% of close 
pairs. Similarly, examination of regions around 
bright stars shows essentially no spurious objects 
in the magnitude range we use. Finally, the over- 
lap between scanlines (roughly 1') means that, for 
galaxies smaller than 20" (which includes essen- 
tially all galaxies in our sample), we can reject ob- 
jects which overlap the edges of the CCDs without 
leaving any gaps. 

We turn now to describe^ our criteria for select- 
ing isolated galaxy pairs. The galaxy angular cor- 
relation function, uj(0), quantifies the ratio of clus- 
tered to unclustered systems; we thus work at sep- 
arations where w(0) > 1, to ensure that the excess 
of pairs over a random distribution is more than a 
factor of two (see Infante, de Mello, & Menanteau 
1996 and Carlberg, Pritchet, & Infante 1994 for 
further discussion). 

The SDSS galaxy redshift survey (cf. Zehavi et 
al. 2001, Stoughton et al. 2001) extends only to 
r* = 17.77, so we use the luminosity function of 
the CN0C2 redshift survey by Lin et al. (1999) 
to determine the expected redshift distribution of 
the sample. We carried out transformations of 
the photometric bands, following Fukugita et al. 
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(1996), for each galaxy type (sec Dodclson et al. 
2001 for further discussion of this point). The pre- 
dicted redshift distribution of galaxies, dN/dz, in 
the range 18 < r* < 20 is plotted in Figure 3, 
and of course is consistent with the distributions 
shown by Dodelson et al. (2001). The predicted 
median redshift is (z) — 0.22 with a dispersion of 
0.09. In what follows, we explicitly assume that 
the luminosity function of galaxies in pairs is the 
same as that of isolated galaxies. 

We select groups of galaxies as follows; below, 
we will limit our analysis to those groups which 
have exactly two members. We examine all galax- 
ies in the catalog, and find all galaxies within an 
angular separation of ^ = 15" (in practice, we also 
require 9 > 2", as the deblender rarely separates 
objects reliably at closer separation). 15" corre- 
sponds to 37 kpc at the median redshift of 
the sample. We then merge all groups which have 
members in common. For each group, we define 
the group radius, Rg, as the angular radius of the 
smallest circle containing the centers of the group 
members (for pairs, Kq is half the separation of 
the two galaxies), and Rjv, as the distance from 
the group centroid to the next galaxy not in the 
group. Our sample of isolated groups is defined by 
the criterion Rtv/Rg > 3, which is based on those 
used for local samples (Karachcntscv 1972: Hick- 
son 1982; Prandoni et al. 1994), who show that 
this decreases the fraction of chance projections 
in group-finding. We note however that our iso- 
lation criterion does not consider galaxies fainter 
than the sample limit of r* =20; in contrast, the 
low-redshift samples referenced above considered 
all galaxies up to three magnitudes fainter than 
the brightest member. We will show below that 
this makes little difference in our analysis. 

This sample contains 15,492 close pairs, but 
only 1175 groups with three or more members, 
too small a sample to measure a robust correla- 
tion function. We therefore carry out the cluster- 
ing analysis only on the pairs. We inspected the 
images of a large number of these pairs. Fewer 
than 3% are spurious detections, e.g. improperly 
deblended bright extended objects or bright star 
spikes, meteor trails, etc.. Fig. 4 shows some ex- 
amples of galaxy pairs from our sample. 

In Table 1 we present the main characteristics 
of the galaxy and pair catalogs. We explicitly as- 
sume that the luminosity function, and thus the 



redshift distribution of galaxies in pairs is the same 
as that of field galaxies. We list the mean redshift 
and eff'ective width of the dN/dz distribution. The 
effective volume, weighted by dN/dz, is listed for 
the samples (equation 6). The correlation statis- 
tics are described below. 

3. The Angulcir Correlation Function of 
Galaxies and Pairs 

We determine the angular correlation function 
of galaxies and of pairs of galaxies using a cata- 
log of randomly distributed points over the survey 
area, and using the estimator suggested by Landy 
& Szalay (1993): 



Ndd -2Ndr + Nr^ 

Nrr 



(1) 



Here Nm is the number of galaxy pairs in a given 
range of separations summed over all SDSS fields, 
Ndr is the number of galaxy-random pairs, and 
Nrr is the number of random pairs. Pair counts 
were done in each scanline separately. The random 
pair counts were obtained from 100 random cat- 
alogs, with galaxies randomly distributed within 
the net survey area. Random counts were normal- 
ized to each scanline separately. We also calcu- 
lated the correlation function using the method of 
Infante (1994); the results arc indistinguishable. 

The resulting angular correlation functions of 
galaxies and galaxy pairs are shown in Figure 5 
and Table 1. The error bars result from 100 boot- 
strap re-sampling variations; these conservative 
errors (see the discussion in Fisher et al. 1994) 
are somewhat larger than the scatter from one 
scanline to another, presented for the galaxy pair 
sample in Figure 6. The galaxy angular correla- 
tion function follows a power law, u> = Au9''^~'^\ 
with 7 = 1.77 ± 0.04, and an amplitude essen- 
tially identical with that given by Connolly et al. 
(2001), when their results for the 18 < r* < 19 
and 19 < r* < 20 slices are combined. In agree- 
ment with the latter paper, we find no breaks in 
the correlation function on scales smaller than one 
degree. 

The galaxy pair correlation function shows the 
same power law behavior as the galaxies, but with 
a significantly larger amplitude. 

The error bars determined by the bootstrap re- 
sampling method in the galaxy correlations are 
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smaller than the symbols. Scranton et al. (2001) 
discuss the angular correlation and the related un- 
certainties in greater detail; our results are consis- 
tent with theirs. We have also carried out a jack- 
knife estimation of the errors, and find consistent 
results with those presented here. 

We fit each of the correlation functions to a 
power law; the results are given in Table 1. These 
fits use the bootstrap errors, but unlike Connolly 
et al. (2001) do not take the covariance into ac- 
count; we may therefore be underestimating our 
error bars on ro somewhat. 

Our galaxy pair sample is defined without red- 
shifts, over a fairly narrow range of apparent 
magnitude. This could give rise to two types 
of systematic error in the estimated clustering 
strength. First, distant clusters of galaxies might 
have only two members bright enough to enter 
the galaxy sample, and therefore could masquer- 
ade as a galaxy pair; as clusters exhibit apprecia- 
bly stronger clustering than do groups (e.g., Bah- 
call 1988), this could bias the results. However, 
at z 0.22, the density of galaxy pairs is much 
higher than the density of clusters. Thus, any 
contamination from background clusters should be 
small. To check this, we searched for galaxies as 
faint as r* =21.5 around each pair, and found no 
significant density enhancements. Thus there is 
no evidence that many of our pairs are the "tips 
of the iceberg" of distant clusters. 

Second, inclusion of pairs in projection will tend 
to systematically decrease the clustering ampli- 
tude. We have carried out simulations of the ef- 
fects of random projection on the clustering of our 
pair catalog. We find that of order / = 30% of our 
pairs are likely to be chance projections along the 
line of sight. The angular correlation function is 
depressed by a factor (1 — /)^ by projection; this 
suggests that the correlation amplitude of groups 
listed in Table 1 should be multiplied by a factor 
of order two. 

4. The Three-Dimensional Correlation 
Function of Galaxies and Pairs 

We relate our lo{9) measurements to the three- 
dimensional correlation functions ^(r) via Lim- 
ber's Equation (Limber 1953). The spatial cor- 
relation function is assumed to be a power law 
weighted by the standard phenomenological evo- 



lutionary factor, (^{r,z) = (^-^^ ^(1 + z)-^^+''\ 
where r is the proper distance, ro is the proper 
correlation length, and e is the clustering evolution 
index, e = 7 — 3 corresponds to clustering fixed 
in comoving coordinates, while e = represents 
stable clustering in physical coordinates (Phillips 
et al. 1978). 

Parameterizing the angular clustering as u){6) = 

^(1-7), 

ro is given by 



^ C gizKdN/dzfdz 



' ~ ^OJ 



AZ' C 



[J^^{dN/dz)dz] 



(2) 



(e.g. Efstathiou et al. 1991; see also Peebles 1980 
and Phillips et al. 1978). C is a constant involving 
purely numerical factors. 



C = 7t'/' 



1/2 r[(7-i)/2] 



r(7/2) 



(3) 



and the function g{z) depends only on e, 7, and 
cosmology: 



g(z) =(^] X^-^Fix)il + z)-(3+^-T) 

\dx J 



(4) 



where x{z) is the coordinate distance and F{x) is 
F{x) = [1 - {Hoaox/cf{Qo - !)]'/'• (5) 

We use a cosmology with qo = 0.1 and A = 0, 
and an evolution index of e = as suggested by 
observations (Carlberg et al. 2000). 

The Ho dependence of g{z) is canceled by that 
in the measured value of ro. The strong depen- 
dence of on dN/dz is clear in this equation. 
Note in particular that the above equations are 
not affected by galaxy evolution, except in the cal- 
culation of the redshift distribution dN/dz (Pee- 
bles 1980, eqs. [56.7] and [56.13]). However, Lim- 
ber's equation, as presented here, does assume 
that galaxy clustering is independent of luminos- 
ity; the high-redshift objects in our narrow mag- 
nitude slice are of higher luminosity than the low- 
redshift objects. Clustering is in fact weakly de- 
pendent on luminosity (e.g., Norberg et al. 2001 
and references therein), but the quantification of 
this effect on our results is beyond the scope of 
this paper. 

Because our sample covers a narrow magnitude 
range, and the galaxies in pairs have essentially the 
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same mean magnitude as the entire galaxy sam- 
ple (Table 1), we expect the dN/dz of pairs to be 
identical to that of galaxies (Figure 3). This al- 
lows us to estimate the correlation lengths, ro, for 
the clustering of galaxies and pairs. We find that 
pairs have a significantly larger correlation length 
than do galaxies; rg = 4.2±0.4 Mpc for galax- 
ies and 7.8 ± 0.1h~^ Mpc for pairs. These values 
are to be compared with ~ 5 Mpc for galax- 
ies (Davis & Peebles 1983; Loveday et al. 1996, In- 
fante & Pritchet 1995), ~ 7-8 h''^ Mpc for groups 
(Carlberg et al. 2001, Girardi et al. 2000, Merchan 
et al. 2000) and ~ 15-20 ft," ^ Mpc for rich clusters 
(Bahcall & West 1992). The errors in our derived 
ro include both statistical uncertainty, as well as 
an uncertainty added in quadrature due to our im- 
perfect knowledge of the redshift distribution (cf., 
the discussion by Dodelson et al. 2001). 

In the following section, we will relate the clus- 
tering length to the mean inter-system separation 
of galaxies and groups. We determine the latter 
from the effective volume of the sample, based on 
the model redshift distribution: 

J^2l^{dN/dz)dz ^ ' 

which then yields the number density, n = 
Nsy^ems ^ mean separation, d = (^)^^"^- We 
find the mean separation of galaxies and galaxy 
pairs is 3.7 and 10.2 h^^ Mpc, respectively. Note 
that random superpositions will artificially depress 
the pairs value of by a factor (1 — /)^/^ 0.9; 
thus the true value of d for pairs is probably closer 
to ll/i-i Mpc. 

5. The ro — d Relation 

The correlation lengths of galaxies and galax;y 
pairs are compared with the global richness- 
dependent correlations of galaxy systems (Bahcall 
& Soncira 1983, Szalay & Schramm 1985, Bah- 
call 1988) in Fig. 7. Here the correlation scales 
ro of various systems are plotted as a function of 
the mean separation of the objects, d; since richer 
systems are rarer, the quantity d scales with the 
richness or mass of the systems. The data points 
include groups and clusters of galaxies of differ- 
ent richnesses and different samples: rich Abell 
clusters (Bahcall & Soneira 1983, Peacock & West 
1992, Postman et al. 1992, Lee & Park 2000; note 



that Richness clusters cannot be included in 
any statistical analyses since they are no complete 
samples of these objects); APM clusters (Croft et 
al. 1997 (C97); Lee & Park 2000 (LPOO); the latter 
find considerably larger ro in their re-analysis of 
the APM clusters); EDCC clusters (Nichol et al. 
1992): X-ray selected clusters (REFLEX survey, 
Bohringer et al. 2001; XBACS survey, Abadi et 
al. 1998, Lee & Park 2000); and groups of galaxies 
(Merchan et al. 2000; Girardi et al. 2000). The 
well-known dependence of ro on d (Bahcall 1988, 
Bahcall & West 1992) is clearly observed. The 
new result presented for pairs of galaxies, with 
d = 10.2 /i-i Mpc and ro = 7.8 ± 0.7 h''^ Mpc, 
fits well within this universal clustering trend; it 
is significantly larger than the galaxy correlation 
scale, and somewhat smaller than the correlation 
scales observed for small groups of galaxies. 

Also presented in Fig. 7 is the expected vq ~ d 
relation for two cosmological models: LCDM 
(n„ = 0.3, nL ^ 0.7, h = 0.7) and Standard 
CDM (SCDM; 17„ = 1, /i = 0.5) (both models are 
normalized to the present-day cluster abundance). 
These relations are obtained from large scale, high 
resolution cosmological simulations (Governato et 
al. 2000; Colberg et al. 2000; Bahcall et al. 2001). 
Figure 7 highlights the fact that the SCDM model 
is highly inconsistent with the clustering data - 
not only for rich clusters, as was previously known 
(e.g., Bahcall & Cen 1992, Croft et al. 1997, Gov- 
ernato et al. 2000, Colberg et al. 2000) - but also 
for small groups and galaxy pairs; the correlation 
strength observed for all these systems is system- 
atically stronger than predicted by SCDM models. 
The LCDM model provides a considerably better 
match to the data; the lower the value of flm, the 
stronger the group correlations predicted. This 
'canonical' model is consistent with most other 
observations of large scale structure, clusters of 
galaxies, supcrnovae la, and the cosmic microwave 
fluctuations {e.g. , Bahcall et al. 1999, Wang et al. 
2000). In fact, the data may even suggest a some- 
what lower Cl„i than 0.3, as some of the best data 
points exhibit stronger correlations than expected 
for 0„ = 0.3. 

Note that the lower rg values for the two rich- 
est (but small) sub-samples of APM clusters are 
due to their exceptionally steep correlation func- 
tion slope: 7 ~ — 3 instead of ~ —2; all other 
sample have a slope of ~ —2. This steepness re- 
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duces To, bringing it lower than the other relevant 
data. 

6. Conclusions 

We use a sample of 330,041 galaxies within a 
stripe of area 278 deg^ centered on the Celestial 
Equator, with magnitudes 18 < r* < 20, obtained 
from SDSS commissioning imaging data. Wc use 
these data to select isolated pairs of galaxies. We 
determine the angular correlation function of the 
galaxies and of the galaxy pairs. We find the fol- 
lowing results: 

• Pairs of galaxies are more strongly clustered 
than are single galaxies. The angular corre- 
lation amplitude of galaxy pairs is 2.9 ± 0.4 
times larger than that of galaxies. 

• The power-law slopes of the two correlation 
functions are the same, corresponding to 7 = 
1.77 ±0.04. 

• We measure lv{0) up to 1 deg scales, cor- 
responding to ~ 9 Mpc at the mean 
redshift of 0.22. No breaks are detected in 
either correlation function. 

• Assuming a redshift distribution from the 
CN0C2 survey luminosity function, we in- 
vert the angular correlations and determine 
a three-dimensional correlation length, vq, 
for each sample. We find ro = 4.2 ± 0.4 

Mpc for galaxies and 7.8 ± 0.7/i~^ Mpc 
for pairs; the latter may be biased downward 
somewhat by projection effects. 

• The mean separation between systems is d = 
3.7 and 10.2 Mpc for galaxies and pairs 
respectively. The correlation lengths fit the 
global ro — d relation observed for galaxy 
systems well (Bahcall & West 1994). 

This work suggests a number of follow-up stud- 
ies. We have defined pairs of galaxies from their 
positions on the sky, and we need to quantify what 
fraction of these objects are pairs at the same red- 
shift. We are in the process of carrying out a 
redshift survey of a subset of the pairs sample, 
which will also be useful in tying down the dN/ dz 
relation. We can also use photometric rcdshifts 
from the multi-band photometry of the SDSS to 



define a cleaner sample of galaxy pairs. Finally, wc 
have used just under 300 deg^ of SDSS data in this 
study; the survey has now imaged over five times 
this area. Thus we are in the process of defining 
samples of richer (and thus rarer) systems from 
this larger area, and measuring their correlations, 
to investigate the ro — d relation between for richer 
groups. 
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Fig. 1. — r* differential galaxy number counts. 
The data are the number of galaxies per magni- 
tude per deg^. The thin solid lines represent the 
counts of 12 individual scanlines. The solid cir- 
cles are the overall number counts in the total 
278.13 deg^. The bold solid line is a power law 
fit to the data. The line can be represented as 
log(A''magdeg^) = 0.46r*— 6.03. The galaxy sam- 
ple is complete to r* = 20.5. 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v5.0. 
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Table 1 

Properties of Galaxies and Pairs (18 < r* < 20) 



Property Galaxies Pairs 

Mean rcdshift (z) 0.22 ±0.1 0.22 ±0.1 

Effective Volume [h'^ Mpc^] 17.8 x 10*^ 17.8 x 10" 

(r*) [SDSS band] 19.34 19.30 

Number 330,041 15,492 

(7) 4.94 ±0.02 (1.77) 13.54 ± 0.07 (1.76) 

Correlation length ro [/i"^ Mpc] 4.2 ± 0.4 7.8 ± 0.7 

Mean spacing d Mpc] 3.7 10.2 



Note.— i/o = 100ft km s-^ Mpc'S qo = 0.1, A = 0, e = 0; 
Area covered: 278.13 deg^ 
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Fig. 2. — Distribution on the sky of galaxies from 
runs 752 and 756 with 18 < r* < 20. The region 
shown is centered on the right ascension range over 
which the two runs overlap. 
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Fig. 3. — Normalized redshift distribution of 
galaxies in the range 18 < r* < 20 based on lu- 
minosity functions from the CN0C2 redshift sur- 
vey (Lin et al. 1999). The median redshift is 
{z) = 0.22, with a FWHM of 0.09. 
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Fig. 4. — Examples of galaxy pairs in r* . The 
cross-hair marks one of the members. Each panel 
is 1.34 X 1.34 arcmin^. These images include fore- 
ground stars as well, and are recreated from the 
SDSS "Atlas Images", postage stamps extracted 
around each detected object. Thus there is no 
background sky noise beyond the boundary of each 
object. 
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Fig. 5. — Angular correlation function of galaxies 
and galaxy pairs for 18 < r* < 20. The x-axis 
is separation in arcsec, measured on a logarithmic 
scale on the lower axis, and a linear scale on the 
upper axis. The open points are the measured 
galaxy correlation function, and the solid points 
are the measured pair correlation function. The 
thin solid line is the fit to the Connolly et al. (2001) 
correlation function with 7 — 1.73, scaled to our 
magnitude range. The dashed line is the best fit 
model for galaxies, with ro = 4.2 Mpc. The 
dot-dashed line is the best fit model for pairs with 
ro = 7.8 Mpc. The error bars are determined 
by bootstrap resampling, and are smaller than the 
symbol size for the galaxies. 
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Fig. 6. — Angular correlation function of galaxy 
pairs for each scanlinc separately. The solid line 
is the fit to the SDSS galaxy correlation function 
from Connolly et al. (2001), scaled to our magni- 
tude range. 
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Fig. 7. — Correlation length tq versus mean sepa- 
ration d (= n~^/^) for galaxies and pairs, as well as 
groups and clusters of galaxies (see § 5, which de- 
scribes the different samples plotted) . The predic- 
tions from simulations of two cosmological mod- 
els, LCDM and SCDM, are shown for compari- 
son, as well as the original approximate relation 
of ro « 0.4 d for rich systems (Bahcall 1988). 
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